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ABSTRACT 

Intergalactic Mg II absorbers are known to have structures down to scales 
~ 10 25 pc, and there are now indications that they may be fragmented on scales 



< 10" 2 - 5 pc (Hao et al., astro-ph/0612409). When a lensed quasar is microlensed, 



the micro-images of the quasar experience creation, destruction, distortion, and 
drastic astrometric changes during caustic-crossing. I show that quasar mi- 
crolensing can effectively probe Mg II and other absorption "cloudlets" with sizes 
~ 10 — 10~ 20 pc by inducing significant spectral variability on the timescales 
of months to years. With numerical simulations, I demonstrate the feasibility of 
applying this method to Q2237+0305, and I show that high-resolution spectra of 
this quasar in the near future would provide a clear test of the existence of such 
metal-line absorption "cloudlets" along the quasar sight line. 

Subject headings: gravitational lensing — intergalactic medium — quasars: ab- 
sorption lines 



Introduction 



Mg II absorbers towards the quasar sight lin es have been syste matically studied since 

Lanzetta et al.l (119871 ) (for more recent studies, see lZibetti et al.ll2005l and references therein) . 

Simila r Mg II absorbers were subsequently seen in gamma-ray burst (GRB) spectra. iProchter et al 
( 120061 ) compared GRB and quasar sight lines, and found a significantly higher incidence to- 
ward the former. They proposed three possible effects to explain this discrepancy: (1) faint 
quasars are obscured by dust associated with the absorbers; (2) Mg II absorbers are intrinsic 
to GRBs; (3) gravitational lensing of the GRB by the absorbers. However, they concluded 
that none of these effects provide a satisfactory explanation. 
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Frank et al.l (120061 ) proposed a simple geometric solution to the puzzle. They argued 
that if Mg II absorption systems are fragmented on scales < 10 16 cm, similar to the beam sizes 
of GRBs, then the observed difference in incidence of Mg II absorbers would simply reflect 
the difference in the average beam sizes between GRBs and quasars, with quasars being 
on average several times as big. This explanation predicts that absorption features due to 
intervening Mg II cloud fr agments should e volve as the size of GRB afterglow changes, which 
has now been observed by lHao et al.l (120061 ) . However, structures of the Mg II ab sorbers down 
to the size of ~ 10 16 cm cannot be directly inferred from their spectral features. iRauch et al. 



(120021 ) put the strongest upper limits on Mg II absorber s izes to date. They observed the 
spectra of three images of Q2237+0305 (IHuchra et al.lll985l ) and found that each line of sight 
contained individual Mg II absorbers at approximately the same redshift, but with distinct 
spectral features. Thus these absorbers are part of a complex that extends at least ~ 500pc, 
but the sizes of the individual "cloudlets" must be smaller than 200 — 300pc based on the 
separation of the macro-images. 

Some, if not all, strongly lensed quasars ar e also Kravitational ly microlensed by the com- 
pact stellar- mass objects in the lensing galaxy (1Wambsganssll2006l ) . and Q223 7+0305 was the 



first l ensed quasar to be fo und to exhibit significant microlensing variability (ICorrigan et al. 
199ll ; IWozniak et al.l |2000| ) . The macro- image of a microlensed quasar is split into many 



micro-images, and when the source moves over the caustic networks induced by the mi- 
crolenses, those micro-images will expand, shrink, appe ar, disappear and experience drastic 
astrometric shifts over timescales of months or years (ITreyer fc Wambsganssl 12004 ) . The 
angular sizes of major micro-images are usually of the same order as those of the quasars, 
and during the shape and position changes of these images, absorption structures of similar 
scale along their sight lines will likely imprint significant variations on the spectrum. 



Brewer fc Lewis! (120051 ) pioneered the theoretical investigation of quasar microlensing as 



a probe of the sub-parsec structure of intergalactic absorption systems. They concluded that 
variation in the strength of the absorption lines over timescales of years or decades caused 
by microlensing can be used to probe the structures of Lyman a clouds and associated 
metal-line absorption systems on scales < O.lpc. However, as I will show, they significantly 
underestimated the relevant timescales for spectral variability given the sizes of the systems 
they considered. Thus, they substantially overestimated the scales of absorption structures 
that microlensing can effectively probe. 

In the following section, I lay out the basic theoretical framework of the method. Then 
in § El I present a numerical simulation of the microlensing of Q2237+0305. I show that 
micro-images of this quasar can be used to probe structures of Mg II and other metal-line 
absorption clouds on scales of ~ 10 14 — 10 16 cm by monitoring the spectral variations of 
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absorption lines over months or years. Finally in § HI I summarize the results and discuss 
their implications. 



2. Varying Microlensed Quasar Image as A "Ruler" 

I begin with a brief summary of notation. Subscripts "/", "s", "o" and "c" refer to 
the lens, source, observer and absorption cloud plane, respectively. The superscript "ray" is 
used to refer to the light ray on the cloud plane to distinguish it from the cloud. The angular 
diameter distance between object x and y is denoted D xy and is always positive regardless of 
which is closer; in particular, D x refers to the angular diameter distance between the observer 
and object x. The vector angular position of object x is denoted 6 X , while its redshift is 
denoted z x . 

Consider an absorbing cloud that is confronted with a "bundle of light rays" making 
their way from the source to the lens to the observer. Let r r c ay be the position vector of 
a point on the plane of the cloud. The line depth (Ax) of an absorption line centered at 
wavelength A is given by: 



(A x )= / ax(rr)A x (rr)d 2 r™v , 



ax(r r c ay )d 2 r r c ay 



where ax{r r c ay ) is the surface density of the "ray bundles" on the plane of the absorption 
cloud with the rays weighted by the surface brightness p rofile of the quasar an d Ax(r™ y ) is 
the absorption fraction at r r c ay for light of wavelength A f lBrewer &: Lewisll2005l ). 



2.1. Basic Geometric Configurations and Motions 

The absorption cloud can be located either between the lens and observer or between 
the lens and source. In the former case, the angular position of the light ray at the cloud 
plane is T c Ay = 0j, so r^ ay = 0{D C . The projected light rays on the cloud plane maintain 
the exact shapes of the quasar images, and their physical extents are proportional to the 
distance to the observer. 

If the cloud is betwee n the lens and the source, it can be easily shown that 6 r c ay is given 



by (IBrewer fc Lewis! 120051 ): 
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so 



r ray 



U\ s J L>i s 



(3) 



The lens-source rel ative proper motion , with time as measured by the observer, is given (in 
other notation) by lKayser et al.l (119861 ). 



1 



Vis 



1 



1 + z s D s 1 + ziDi l + zi DiD s 



(4) 



where v a , vj, v are the transverse velocities of the source, lens and observer, relative to the 
cosmic microwave background (CMB). In particular, 



VCMB — (VcMB ' Z J Z ! 



(5) 



where z is the unit vector in th e direction of the lens and vcmb is the heliocentric CMB 
dipole velocity (IKochanekl 120041 ) . 



The formulae in this subject can be greatly simplified by proper choice of notation. To 
this end, I defined the "absolute" proper motion of an object x moving at transverse velocity 
v x to be: 

1 + z x D x 

I also define the "reflex proper motion" of an object x relative to the observer-lens axis to 
be: 

, 1 v Di x 



sgn^a, - z^ 



zi D t D a 



(7) 



Then equation (jlj) can then be simplified, 



Vis — Vabs,s Mabs,/ + V 0) l,s- 

Similarly, the lens-cloud relative proper motion is given by: 

Vic = Mabs.c — Vabs.l + Vo.l.c 



(9) 



Note that the last term has a different sign depending on whether the cloud is farther or 
closer than the lens (see eq. [7j). 



2.2. Bulk Motion of the Un-microlensed "Ray Bundles" 

Equations (j2j), (jSJ) and (jUJ) are the key formula needed to carry out the simulation in 
|3j In this section, I discuss the bulk motion and relevant timescales of the macro-image 
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and its associated "ray bundles" on the cloud plane for the underlying case that the image 
is not perturbed by microlenses. 

When the source moves at fj, ls , the angular positions of the rays that compose the images 
are also changing with respect to the observer-lens axis. If the macr o-image is unperturbe d 



by the microlenses, the relative proper motion fj, H is simply given by (IKochanek et al.lll996l ). 



Hu = M- » ls , (10) 

where Ai is the magnification tensor. At the same time, the intersection of "ray bundles" 
with the cloud plane also change their angular positions relative to the lens. If the cloud 
plane is between the lens and the source, then from equation (j2J), the "ray bundle"-lens 
relative proper motion /x™ y is a linear combination of fx ls and /% weighted by distances, 

Substituting equation ( flOl) into equation ( flTi) yields the bulk proper motion of the un- 
microlensed "ray bundle" relative to the lens, 



./ay 



_ D lc D 8 \ D lc D a _ 
D ls D c J + D ls D c 



Vis- (12) 



where X is the unit tensor. 



When the "ray bundles" are between the source and the lens, their relative bulk proper 
motion is simply, 

»Zuik = »u = Mfr) ■ (13) 

There are two major effects that may induce variability of absorption lines toward a 
lensed quasar. One is the creation, destruction, distortion and astrometric shifts of micro- 
images, which probe structures similar to the size of the micro-images. The other is the bulk 
motion of the "ray bundles" relative to the cloud. On the cloud plane, this motion has an 
angular speed of A/x Zc bulk = ^ y bulk — £t ic , and the time t cc required for the "ray bundles" to 
cross a cloud of transverse size R c is, 

^ cc = n I a \ ~ n \' (1^/ 

JJ c\ L ±H'lmlk,lc\ LJ c\Hbulklc~ A*/cl 



3. Application to Q2237+0305 



In October 1998, iRauch et al.l (120021 ) obtained high-resolution Keck spectra of images 
A, B and C of Q2237+0305. They found Mg II absorption lines at redshifts of z — 0.5656 
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and z = 0.827 in the spectra of all three images, but absorption profiles of the individual 
sight lines differed (e.g., Figs. 6 and 10 of their paper). Therefore, they concluded that the 
Mg II complexes giving rise to these absorption features must be larger than ~ 0.5kpc, while 
the individual Mg II components must be smaller than ~ 200 — SOOfc^pc. Q2237+0305 is 
also one of the most observed and studied lensed quasars with obvious microlensing features, 
and the properties of the system are well known. These factors make it an ideal object to 
investigate. 



The comprehensive statistical study of this lens by IKochanekl (120041) showed that the 



Mortonson et al 



(b()05h demonstrated that 



size of the quasar is ~ 10 15 /i _1 cm — 10 16 /i _1 cm. 
the source size has a much more significant effect on microlensing models than the source 
brightness profile. For simplicity, in the simulation, I model the source as uniform disks 
with four different sizes: 10 15 /i _1 cm, 3 x 10 15 /i _1 cm, 5 x 10 15 /i~ 1 cm and 10 16 /i _1 cm. Uniform 



grids of rays are traced from image plane to source plane (jWambsganssI Il990l ) . Because 



structures of interest have similar sizes as the source, finite-source effects must be taken into 



accou nt. The grid size used has an angular scale 1/10 of the smallest source. iKochanek 
( 120041 ) demonstrated that a Salpeter mass function cannot be distinguished from a uniform 
mass distribution and found the mean stellar mass to be (M) ~ O.O37/i~ 1 M . For simplicity, 
I assign all stars in the simulation the same mass of O.O4/i -1 M . Rays are shot from a region 
extending 47(#e) on each side, where (#e) is the Einstein radius of a .04/t~ 1 M f 7) star. I 
adopt a convergence and shear for image A of (k, 7) = (0.394, 0.395) from IKochanekl (120041 ) . 
and set the stellar surface density Four different trajectories, oriented at 0, 30, 60 

and 90 degrees with respect to the direction of the shear are studied. Positions on both the 
image plane and the source plane are recorded once a ray falls within a distance of 2 times 
the largest source size from any trajectory on the source plane. A total length of 5(#e) along 
each trajectory is considered. 

Throughout the paper, I adopt a ACDM cosmology with Q m = 0.3, = 0.7 and 
H = 100 /tkms" 1 Mpc" 1 . The lens and source are at redshifts z\ = 0.0394 and z s = 
1.695 Jrluchra et alJ ll985h. These imply (D s , D h D u ) = (1223, 113, 1180) h~ x Mpc. Based on 
Kochanek J2004h . I adopt transverse velocities of the lens, source and observer of (vj, v s , v ) = 
(300, 140, 62)kms~ 1 . The lens and source absolute proper motions and the source reflex 
proper motion are (// a b s ,^ A*abs,s > Ho,i,s) — (0.54, 0.009, 0.11) h /zasyr -1 . So the lens absolute 
proper motion dominates the lens-source relative proper motion. It can easily be shown that 
the absolute and reflex proper motions of the cloud are much smaller than the absolute proper 
motion of the lens unless the cloud redshift is close to or smaller than the lens redshift. For 
Q2237+0305, the lens redshift is very small compared to the source, so most likely z c ^> z\. 
Thus in following analysis, I focus on cloud with redshift z c > zi, and ignore absolute and 
reflex proper motions of the source and the cloud. In this case, the source and the cloud 
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share the same relative proper motion, 

Vic = Vis = -A*abB,i ( 15 ) 

In the simulation, as a practical matter, I hold the positions of the observer, lens galaxy 
(as well as its microlensing star field) fixed, and allow the source to move through the source 
plane at fi ls . Then by equation (j!5J) . the cloud has the same relative proper motion as the 
source. At any given time, the angular positions of the "ray bundles" are calculated using 
equation (j2J). Then by simply subtracting the angular position of the source at that time, 
the "ray bundles" positions are transformed to the reference frame of the cloud. 



During short timescales, micro lensing causes centroid shifts of the macro- image (ILewis fc Ibata 



19981 ; iTreyer fc Wambsganssll2004j ). with respect to the steady bulk motion of "ray bundles" 
relative to the cloud, which is described in § 12.21 If the direction of relative lens-source 
proper motion is the same as the lens shear, then by subtracting [i ls from equation ffl2|) . one 
finds that the "ray bundles" have their maximum bulk proper motion relative to the cloud, 



\ h AcAA ( 1 A 

D ls D c J \1 - K-7 / 

1 - 1.037^J hfiasyi' 1 . (16) 

If the source moves perpendicular relative to the lens shear, A/i Zc = [1— D LC Ds/(D L sDc)]\l/(l- 
k + 7) — l\^is, which is approximately for the (k, 7) of image A. 

Figure [JJ and Figure [2] show the results for a source trajectory that is parallel to the 
shear direction. The bottom panel of Figure [1] shows the magnification pattern on the source 
plane, with a series of 4 concentric circles centered at 5 source positions; and the top panels 
show the images at these positions relative to the source (which has the same proper motion 
as the cloud). Different colors represent different source sizes. The middle panel of Figured] 
shows the light curves for the 4 source sizes with the blue dash lines used to mark the times 
for the five source positions. Figure [2] shows the "ray bundles" positions in the cloud frame 
at redshifts 1.69, 0.83, 0.57, 0.1. The 5 different columns show the 5 positions corresponding 
to those in Figure [JJ 

In the top row of Figure [2J one can see that the "ray bundles" at z c = 1.69, which is 
very close to the quasar, have almost exactly the same size and shape as the source and 
that the bundles show almost no bulk motion. The density of "ray bundles" clearly have 
the imprints from magnification pattern shown in the bottom panel of Figure [TJ So if an 
absorption cloudlet has a similar or somewhat smaller size than a source that is sitting 



- 8- 



directly behind it, the depth of its corresponding absorption line will change dramatically as 
the source crosses the caustics. The magnification close to a fold caustic is proportional to 
the inverse square root of the distance from it, so for a cloud with angular size 9 C < 6 S , the 
fractional change in absorption line depth caused by the caustics scales as (# c /#s) 3 ^ 4 - Hence, 
for cloud close to the quasar redshift, structures on scale ~ 10 14 — 10 16 /i -1 cin (depending on 
the source size) will be probed over few-month to few-year timescales (i.e., the timescale of 
typical caustics crossings). 

The fourth row of panels of Figure [2] shows "ray bundles" for z c = 0.1, which is close to 
the lens redshift. A distinct difference between these "ray bundles" from the ones at z c = 1.69 
is that they are split into many groups of bundles, which correspond to the micro-images 
in the top panel of Figure [TJ I dub these groups of "ray bundles" as micro-images in the 
following discussions. Most of these micro-images are stretched one-dimensionally, and most 
rays are concentrated in a few major micro-images. The micro-images in different columns 
have drastically different morphologies and positions. If there are cloudlets of similar sizes 
as these micro-images distributed on the cloud plane, then the absorption spectra will show 
multi-component absorption features at any given time. These components will experience 
drastic changes in line depth, with some components disappearing and other new components 
appearing during the course of months or years, as the source crosses the microlensing 
caustics. Another important characteristic is that the bulk of these micro-images are moving 
in the same direction as the source. This motion is described by equation ( fl6l) . which yields 
~ 0.8/z/zasyr. So from equation (JHJ), structures as large as ~ 3. On x 10 16 h~ x cm will 
be crossed in n x 10/i -1 yr by the bulk motion of the "ray bundles". Hence, the effects 
caused by the micro-images and the bulk motion of the bundles together probe scales of 
~ 10 14 — 10 16 /i~ 1 cm on timescales of months to years. 

The second and third row of panels in Figure [2] refer to clouds at intermediate redshifts 
between the lens and the source. The ir redshifts, z r = 0. 83 and z c = 0.57, are close to the 
Mg II absorption systems observed by iRauch et al.l (120021 ). As expected, the characteristics 
of these "ray bundles" are intermediate to those shown in rows 1 and 4. The overall shapes 
of the micro-images are close to that of the source, with magnification patterns imprinted 
on them. And they also clearly show multiple components, which appear and disappear 
as the source crosses caustics. The angular sizes of the images are close to that of the 
(unmagnified) quasar. These micro-images could probe clouds with angular size from a 
factor of few smaller to a factor of few larger than the source size, which corresponds to 
scales of ~ 10 14 — 10 16 /i _1 cm. From equation (j!4p and (1161) . the bulk motion of the bundles 
will probe cloud structures ~ 0.8nx 10 16 h^ 1 cm and ~ 1.3n x 10 16 h^ 1 cm during fix 10 h^ 1 yr 
for z c = 0.83 and z c = 0.57, respectively. 
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For a source trajectory that is perpendicular to the shear, there will be almost no bulk 
motion relative to the cloud for the "ray bundles" , while the magnitude of bulk motions for 
intermediate trajectories is a fraction of the parallel case depending on their angles relative to 
the shear direction. And the micro-images of these trajectories share similar properties with 
those of the trajectory that is parallel to the shear. Figure [2] shows that micro- image motions 
are in the same orders as the bulk motion of macro-image over the timescales considered. 
Therefore, trajectory direction only has a modest impact on the cloud sizes probed. 



These results are in great contrast with those of iBrewer &: Lewis! ( 120051 ). who claimed 
quasar microlensing for image A of Q2237+0305 can induce considerable variability of ab- 
sorption lines associated with structures as large as O.lpc during the course of years to 
decades. According to their analysis, the effect is largest when the cloud is very close 
to the source, and for example, the timescale of line strength variation for a O.lpc cloud 
very close to the source is given as ~ 16.2yr. However, in their analysis, they effectively 
assumed the relative lens-cloud proper motion //j c = 0. This would lead to a timescale 
t cc ~ R C D[/ iyiD s ) ~ 16(i? c /0.1pc)yr for absorption cloud near the source redshift (they 
adopted v^ = 600km s" 1 ), which is in agreement with column 4 of their Table 1. In fact, I 
showed that, when peculiar velocities of the source and cloud are ignored, fj, lc = n ls (eq. !T5]) . 
which is not negligible. Even considering realistic peculiar motions of the source and the 
cloud, it still leads to time scales that are more than one order of magnitude slower than 
those predicted for fii c = 0. In addition, when the cloud is close to the source, the angular 
sizes of the clouds they considered are orders of magnitudes larger than their source size, so 
effects of changes i n magnification pattern on the source plane alone have very little impact 



as well. Therefore, IBrewer fc Lewis! ( 120051 ) significantly overestimated the cloud size to which 



microlensing is effectively sensitive. 



4. Discussion and Conclusion 

I have shown that there are two effects that might induce variation of absorption lines 
along the sight lines to lensed quasar. One effect is caused by the drastic morphological and 
positional changes of micro-images when the source crosses the caustic network. The other 
effect is due to the bulk motion of the "ray bundles" relative to the absorption clouds. I have 
laid out a basic framework in studying these effects for microlensed quasars in general. And in 
particular, I perform numerical simulations to apply the method to image A of Q2237+0305. 
I demonstrated that the combinations of these two effects probe 10 14 — 10 16 h~ x cm absorp- 
tion cloudlets between the lens and the source over timescales of months to years. The 
existence of these cloudlets will be revealed by either changes in line depths or appear- 
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ances/disappearances of multi- absorption components. Spectra should preferably be taken 
during the course of caustic crossings, which can be inferred fr om photome t ric m onitoring; 
programs of lensed quasars. In fact, the Mg II lines observed by lRauch et al.1 (120021 ) about 8 
years ago already show different multi-components along sight lines of three different macro- 
images, implying they might be caused by fragmented cloudlets with similar sizes as the 
micro-images. A similar high-resolution spectrum taken in the near future would provide 
a definitive test of the existence for structures of Mg II or other metal-line absorbers at 
the scales of 10 1 4 — 10 16 h~ x cm. I f the spectral variations are indeed observed, a statistical 
study similar to iKochanekl (120041 ) will be required to infer the properties of the cloudlets. 
Moreover, a time series of spectra may provide additional constraints to quasar microlensing 
models. 
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Fig. 1. — Caustics network (bottom panel), light curve (middle panel) and images relative 
to the source position (top panel) for a source trajectory that is parallel to the lens shear 
(the direction of x-axis). Different colors represent the 4 different source sizes. 
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Fig. 2. — Physical positions of "ray bundles" in the frame of cloud at redshifts 0.1, 0.57, 
0.83 and 1.69. The 5 different columns correspond to the source positions shown in Figure [TJ 
x-axis is the direction of the lens shear. 



